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1 .  INTRODUCTION 


The  collective  bremsstrahlung  recoil  force  on  a  relativistic  test  particle 
participating  in  a  brensstrahlung  process  in  a  nonequilibrium  relativistic 
bean-plasna  system  is  given  by1'2 

*S  =  -/  <2*)-9  d3£e  d3*  d3jt  t[t  -  t].(»  fe) 

P0  (i ) 

x  +  ^ • 

Here*  v$|  +^($,lt)6(u)jij!  “  +  (^  “  bremsstrahlung  transition 

rate  (probability  per  unit  time)  for  scattering  of  particles  of  species  a  and 

3,  velocities  vQ  and  v^,  and  momenta  pa  and  p^,  where  the  scattering  produces 

momentum  transfer  t  and  the  emission  of  a  photon  in  mode  o  with  wave  vector  & 

and  frequency  w%i  N$  is  the  photon  (plasmon)  distribution  function;  and  f+ 

P$ 

is  the  charged  particle  distribution  function.  Here  and  throughout,  the  units 
are  chosen  such  that  ft  -  1.  In  this  report  equation  (1)  is  derived  from  first 
principles.  This  expression  is  useful  in  obtaining  an  expression  for  the 
collective  bremsst rah lung  transition  rate  by  direct  comparison  with  another 
expression  for  the  collective  bremsst rah lung  recoil  force,  which  is  determined 
from  the  equation  of  motion  for  a  dynamically  polarized  test  particle  under¬ 
going  bremsstrahlung. 1 ' 2  This  transition  rate  is  important  in  calculations  of 
collective  bremsstrahlung  and  the  conditions  for  the  occurrence  of  the  brems¬ 
strahlung  radiative  instability  in  relativistic  beam-plasma  systems  (see  the 
work  of  Akopyan  and  Tsytovich,  Selected  Bibliography). 

In  section  2  the  particle  balance  equations,  including  bremsstrahlung  and 
inverse  bremsstrahlung  processes,  are  derived  in  terms  of  the  basic  brems¬ 
strahlung  transition  rate.  In  section  3,  the  soft  photon  or  quasiclassical 
limit  of  these  equations  is  obtained  and  separated  into  spontaneous  and 

V.  N.  Tsytovich,  Bremsstrahlung  of  a  Relativistic  Plasma,  Tr,  Fiz.  Inst. 
Akad.  Sauk  SSSR,  66  (1973),  191  [Proc  (TRUDY)  of  P.  S.  Lebedev  Physics  Inst., 
66_  ( 1975),  199]. 

2A.  V.  Akopyan  and  V.  N .  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plasmy,  1  (1975),  673  [Sov.  J.  Plasma  Phys.,  1_  (1975),  371], 


induced  contributions.  Expressions  are  obtained  for  the  diffusion  coefficient 
and  induced  recoil  force  due  to  induced  bremsstrahlung.  In  section  4,  the 
general  form  of  the  induced  recoil  force  is  determined  in  terms  of  the  brems¬ 
strahlung  radiation  field,  particle  distributions,  dielectric  permittivity, 
and  bremsstrahlung  transition  rate.  Section  5  summarizes  the  main  results  of 
the  derivation. 


2.  PARTICLE  BALANCE  EQUATIONS 

He  here  derive  the  particle  balance  equations  in  terms  of  the  basic  brems¬ 
strahlung  transition  rate,  for  beam-plasma  configurations  in  which  bremsstrah¬ 
lung  and  inverse  bremsstrahlung  are  the  dominant  processes.  For  this  purpose, 
consider  the  elementary  bremsstrahlung  process  depicted  in  figure  1 .  Here  a 
particle  of  species  a  and  momentum  $a  scatters  off  a  particle  of  species  3  and 
momentum  and  emits  a  photon  of  momentum  &.  Figure  1  is  not  a  Feynman 
diagram  since  it  represents  the  probability  for  the  process  and  not  its  ampli¬ 
tude.  If  the  momentum  decrease  of  particle  a  is  denoted  by  <,  then  the  momen¬ 
tum  of  particle  a  after  the  collision  is  given  by 

-  *  *  (2) 


Then,  by  conservation  of  momentum, 
sion  is  given  by 


momentum  of  particle  3  after  the  colli- 


+  £  -  k  .  (3) 


Figure  1 •  Elementary 
bremsstrahlung  process  which 
removes  particles  of  species 
o  from  point  $a  in  momentum 
space. 
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The  quantity  wpa,$g(^'£)  is  the  probability  per  unit  time  that  in  the  scatter¬ 
ing  of  a  particle  of  the  species  a  and  initial  momentum  £a  off  a  particle  of 
species  0  and  initial  momentum  a  photon  of  momentum  ic  in  mode  o  is  emit¬ 
ted,  and  there  occurs  momentum  transfer  ic  -  k  to  particle  0.  By  time-reversal 
invariance,  the  transition  rate  for  the  inverse  process  (in  which  a  particle 
of  momentum  5a  absorbs  a  photon  of  momentum  ic  and  scatters  off  a  particle  of 
momentum  with  momentum  transfer  i?  to  particle  a)  is  then  given  by 
wpa+K,^g-ic4it^ 11118  Process  is  depicted  in  figure  2  in  the  leftmost 
figure,  which  is  the  time-reversed  process  of  the  figure  to  its  right;  there¬ 
fore,  the  transition  probabilities  of  the  two  processes  are  equal*  The 
inverse  process  is  denoted  by  wfa+£, <,£)  in  the  notation  of  figure  1. 
The  equality  of  the  probabilities  of  the  direct  and  inverse  processes  due  to 
time— reversal  invariance  is  also  known  as  the  principle  of  microscopic  revers- 
ibility  or  reciprocity,  and  leads  to  the  principle  of  detailed  balance*  The 
processes  of  figures  1  and  2  both  deplete  the  density  of  particles  of  type  a 
at  point  in  momentum  space  by  adding  or  subtracting  momentum. 

Figure  2.  Inverse  bremsstrahlung 
process  which  removes  particles  of 
species  a  from  point  pa  in  momentum 
space.  First  equality  follows  from 
time-reversal  invariance. 

■  v'#o+fc|T0-.T+k£® 

The  particle  density  at  point  increases  because  of  the  following  two 
processes.  The  inverse  process  to  that  in  figure  1  produces  gain  at  point 

By  time-reversal  invariance,  this  is  again  given  by  wj?  as 

u  Pa#Ppv 

depicted  in  figure  3.  Similarly,  the  time  reversal  of  the  process  in  figure  2 
also  adds  particles  at  point  $a  with  transition  rate  wpa+£,£g-ic+£(^'^)  <  as 
depicted  in  figure  4. 
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Figure  3.  This  process  adds  particles  at  point 
pQ  in  a  momentum  space. 


Pa  +  * 


P  a  +  * 
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Figure  4.  This  process  also  adds  particles  at 
point  in  momentum  space. 


Taking  these  four  processes  into  account  including  both  spontaneous  and 
induced  emission,  we  find  that  the  particle-balance  equation  giving  the  time 


rate  of  change  of  the  particle  distribution  f>  for  particles  of  type  a  at 

Pot 


point  $n  in  momentum  space  is 


-5 r-'l  <2”-9  d3k  d3?e  d3;  K  *  ') 


-  /  (2»)-9  d*  d^6  d3*  Bkf50,f5BW|a4,56-J<-it^',!) 


♦  /  (2.)-’  d*  dSfj  d3f  »|  £5a.f  «|a,?6(*.*) 


*  /  (2i»-9  as*  d3fs  d3t  fjarf  £j64tj  *  ')  • 


Pa+,c'Pf5" 


The  first  term  in  equation  (4)  represents  the  rate  of  decrease  of  particle 
density  at  point  pQ  due  to  the  bremsstrahlung  process  of  figure  1 .  The  func¬ 
tion  f*  is  the  distribution  function  for  particles  of  species  B,  and  N?  is 

o 

the  photon  distribution  function.  The  factor  N+  +  1  takes  account  of  both 
induced  and  spontaneous  emission.  An  integral  over  all  possible  photon  wave 
vectors,  scattered  particle  momenta  Pg,  and  momentum  transfers  <  also  neces¬ 
sarily  appears  in  equation  (4).  The  factors  of  2*  arise  from  counting  quantum 
states.  Thus,  for  example,  d3k/(2ir)3  ia  the  number  of  quantum  states  per  unit 
volume  with  momentum  4$k  in  the  interval  d3k.  The  same  phase-space  normaliz¬ 
ation  for  the  particle  distribution  function  as  in  Tsytovich1  is  used  here, 
namely. 


n  a  m  J  - r  > 

(  2tt  ) 3  Pa 


where  na  is  the  number  of  particles  of  species  a  per  unit  volume.  The  second 
term  in  equation  (4)  represents  the  rate  of  decrease  due  to  the  inverse  brems¬ 
strahlung  process  of  figure  2.  The  third  term  represents  the  rate  of  increase 

V.  N.  Tsytovich,  Bremsstrahlung  of  a  Relativistic  Plasma,  Tr.  Fiz.  Inst. 
Akad.  Nauk  SSSR,  66_  (1973),  191  (Proc  (TRUDY)  of  P.  N.  Lebedev  Physics  Inst., 
66  (1975),  199). 


due  to  the  inverse  bremsstrahlung  process  of  figure  3.  The  fourth  term  repre¬ 
sents  the  rate  of  increase  due  to  the  bremsstrahlung  process  of  figure  4. 


Combining  terms  in  equation  (4)  produces 


d3k  d3pg  d3< 


[f+  f+  (n 

a  n  a 

it 

f+  -► 

f  ♦  -►  .♦"I 

w* 

g 

[  Pa  P6 

Pa-< 

pg+ic-kj 

r 

(6)  | 

d3k  d3pg  d3<  o  +  r  0  0 

(2ir)9  Pa+</PB-K+k^  [  Pa  Pg  *  Po+lf  Pg-K+k1  k 


Breaking  equation  (6)  into  its  induced  and  spontaneous  parts  results  in 


>f+  /  3f+  \  /  3f+  \ 

Pa  (  Pa)  f  Pa) 

3t  y  3t  Ji  y  3t  J s 


where  the  induced  part  ( 3f +^/3t)^  is  that  proportional  to  the  photon  number 
density  N?  and  the  spontaneous  part  ( 3fpa/9t)g  is  that  independent  of  N+. 
Thus 


GS. 


,  d3^  d3pg  d3lc  0 

-< - 75775  “it 


+ 

Pa*P3 


(«)[< 


f-*-  f+  - 
.  Pa  P0 


f*  * 
Pa~K 


f ♦  + 

pg+K- 


+  WIo^,P6 


LCPafPB  "  fPa+^  fP 


Pg-*+*]} 


(9) 


3.  SOFT  PHOTON  AND  SMALL  MOMENTUM  TRANSFER  APPROXIMATION 

We  next  assume  that  the  momentum  transfer  £  and  the  momentum  (  of  the 
radiated  quantum  are  small  relative  to  the  relativistic  particle  momenta  pa 
and  Pg,  thereby  enabling  Taylor  series  expansions  for  wpa+x,pg-K+jt'  fPa“^' 
fp$+K~i*'  fPa+*'  an<*  equation  (8).  Thus  we  assume  that 

(m«u « iiU'i^n  ■  <>°> 


This  is  evidently  consistent  with  the  Bom  approximation  for  plasma,  namely 
that  the  relativistic  particle  momentum  is  much  greater  than  the  electromag¬ 
netic  impulse  received  by  a  plasma  particle  in  a  time  interval  given  by  the 
inverse  plasma  frequency. 1-5  First  expanding  about  pa  in  the 

first  variable,  we  obtain 


N.  Tsytovich,  Bremsstrahlung  of  a  Relativistic  Plasma,  Tr.  Fiz.  Inst. 
Akad.  Nank  SSSR,  66_  (1973),  191  [Proc  (TRUDY)  of  P.  N.  Lebedev  Physics  Inst., 
§6.(1975),  199). 

^A.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plasmy,  (1975),  673  [Sov.  J.  Plasma  Phys.,  1_  (1975)  371). 

%.  E.  Brandt,  Nonlinear  Dynamic  Polarization  Force  on  a  Relativistic  Test 
Particle  in  a  Nonequilibrium  Beam-Plasma  System,  Harry  Diamond  Laboratories, 
HDL-TR-1994  (1983). 

**H.  E.  Brandt,  Nonlinear  Force  on  an  Unpolarized  Relativistic  Test  Particle 
to  Second-Order  in  the  Total  Field  in  a  Nonequilibrium  Beam-Plasma  System, 
Harry  Diamond  Laboratories,  HDL-TR-1995  (1983). 

*h.  E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi- 
librium  Relativistic  Beam-Plasma  System,  Harry  Diamond  Laboratories,  HDL-TR- 
1996  (1983). 
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(11 ) 


+  k.  .  --  vS  *  *  t 
Pa+^Pp-*-*  Pa»PB-|C+lc  1  3pai  Pa'Pe-*4* 

32W?  ♦  +  + 

1  Pa'PP"*4* 

+  —  <.  K  .  - - - - -  . 

2  i  J  3Pai3Paj 


Next  expanding  equation  (11)  to  second  order  about  Pg  in  the  second  variable 
gives 


W?  ♦  ♦  * 

Pa+K'Pp-*-* 


w£  ■*■  +  (k<  -  <i)  v  —  W* 

Pot'P$  k  1  J  3PSi  Pa»P$ 

i  32wS  -►  3W$  -► 

2  v  1  in  J  3p0i3P3j  1  3Pod 

32W+  32W  +  + 

,  ,  Pa*P3  1  Pa»P3 

+  ^iC  j  Kj)  3pai3pgj  +  2  KiK3  3pai3paj 


Also  in  equation  (8)  the  following  Taylor  series  expansion  applies: 


[Vfe  '  ‘fa-t  V«] 


f>  f+ 
Pa  P0 


3f+  a2f*  "| ' 

fia  -  Kl  *  *  ‘1Ki  [fPs  l13’ 


3fp  1  32fPP 

+  3p^  +  1  ^  3pW3p3n_ 


or  combining  terms  produces 


A  .*•  A.  ' 

•Y  > 

ay  a  .nv  y- 


32f+ 

P0 


'  2  "  *m)f+a  3p$A3pgin  +  Ki  3pJ 


3f+  3f+ 

/  ,  \  Pq  P3 

+  *ci(Kt  1)  3pai  3pgA 


d2f-*- 

Pa 


2  *iKi  apcd^Paj  Pp  * 


if* 

(ki  -  *i)f+  T---- 

1  Pa  3Ppi 

32f+  if* 

*  (*t  -  **)(*.  -  -  <±%  s“ 


it*  it* 

Pa  Pp  1 


r .  \  tra  fp  1  _  fa 

1  1  K£  3Pai  *PPi  2  ICi'C^  Pp  8Poi3Paj 


Substituting  equations  (12),  (14),  and  (15)  in  equation  (8)  produces 


3f+ 

pg 

3p$i 


32f+ 
P5 


3f* 
Pa 


“  2  (K*  “  ki)(*cn  ‘  *m)fp0  3pg*3pgn  +  Ki  3pai  fp0 

3f+  3f+  32f* 

f  u  1  Pa  PS  1  Pa  . 

Ki(KA  i)  3poi  gp^  2  IClKj  3Pai3Paj  P0 


3f* 
P(S 


32f* 


-  (ki  -  Ki)f£a  -  \  (kA  -  KjtXk*  -  KB)f+ 


Pp 


Pa  8P0i8P0m 


-  ‘if£ 


df* 

Pa 


3f+  d£+ 

Pa 


e*  d2f*  1 

(.  \  ^a  P0  1  ,  Pa 

-  <i(k*  -  «i)  3pai  jpw  -  2  **-*11%  3poi3pajJ 

3W$ 


P0  3Pai 


(16) 


3W?  ♦  I"  3£+  3f+ 

♦  -i^r  K  - “>(*.  -  wt  -  c*  - ■*)■**  Pa 


P0  8Pam 


] 


3w^  ♦  r 

Pa>Pp 

+  3Pai  L 


if* 

'itk»  -  k-)fp.  Tth  - 


,  !5.1)  . 

P0  8PaiJj 


The  £irst  and  sixth  teres  of  equation  (16)  cancel.  Also,  the  third  and  eighth 
terns  cancel.  Containing  the  fifth  and  tenth  terns,  the  second  and  seventh 
terns,  and  the  fourth  and  ninth  terms,  and  integrating  the  eleventh  and 
twelfth  terns  by  parts,  we  find  that  equation  (16)  becomes 


(  Pal  t  d3jt  d3Pg  d3,c  a  )  a 

YTT/i  •  •  1 - |wPa'Ps  <,le  [_”  K* '  k‘)(l<*  '  k* 

if*  if*  i2f*  1 

+  Ki  C  jPai  3pgjj  KiKj  apaiSpcj  pgj 

r  32fp 

*  VLit  P  '  '*■ Kk"  '  “*)£?a  tesfe 

if*  if*  “I 

^  Pg  Pa 

+  (ki  '  Ki)<m  3Pgi  3PamJ 

3W$  r  3f  ♦  if*  1) 

Pa'Pg  \M  Pg  ,  Pa 

*  -SST  htk"  '  k")£?«  ^  ’  kik‘£5e  55Sj|  • 


<3~Z* 

)f* 

Pa  3Pgf3Pg 


The  first  and  fourth  terms  cancel;  the  second  and  fifth  terms  combine; 
after  terms  are  rearranged,  equation  (17)  becomes 


/3fPa\  r  d3£  d3Pg  d3£  f  3WPa'Pg 

\“3t"/i  - 1  — tz?r-  M 


3f+ 

O  r  P« 

N+  f*  -r - 

*  Pp  3Paj 


3zf+ 

Pa 


+  K^IC^wS  ♦  (K,k)N>  f+  — 

1  3  Pa'Pg  k  Pg  3Pai3Paj 


-  *i(*j  - 


3f+  3W®  ♦  (ic,k) 

v.'l  -II  Pa>Pgv  J 

j '  3pgj  Spaj. 


x*f* 

*  Pa 


Equivalently, 


/3fpg^  3_  {(",  d3*  d3x  d3pg  *  *  a  ,  ,  o  1 

\  3t  /i  “  apai  |[^  (2*)9  fPBj 


3  ^  ,  d3S  d3S  d3£S  , 

+  ^  ^  - 


if* 

s  \  Pp  0 
3ppj  Pa 


4.  INDUCE}  RECOIL  FORCE  DUE  TO  BREMSSTRAHLUNG 


Equation  (19)  may  be  rewritten  in  the  £orm  of  a  diffusion  equation  as 
follows: 

(-sr)l  -  3poi  (°^3  5i§)  *  Tit  (*.  fJ«)  ■  (2<” 


where  the  induced  diffusion  coefficient  D®,^  is  given  by 


o  d3*  d3x  d3pg  0  r*  *s  a 

Dol3  ■  1  IZW  75TJT  ,CiKiwS„.peC,c'k)sS  % 


and  the  induced  dynamic  friction  force  or  bremsstrahlung  recoil  force  is  given 


Equation  (22)  is  the  sought-after  expression  for  the  recoil  force  on  particle 
a  due  to  the  induced  bremsstrahlung  process*  Energy  conservation  can  be 
explicitly  factored  out  of  the  transition  rate  appearing  in  equation  (22)  as 
follows.  The  bremsstrahlung  transition  probability  "pcpgt*'^)  in  equation 
(22)  must  conserve  energy.  Therefore 


,  VW  WV  V  VW  V 


•* .  *\  *\  •  . 


where  e±  denotes  the  total  energy  of  the  particle  of  species  a  and  momentum 
Pa 

pa.  Using  the  conditions  in  equation  (10)  produces,  to  lowest  order, 

e±  ±  =  e±  -  ic»7_  e±  (24) 

Pa"*  Pa  Pa  Pa 

and 


C?S+t-t  ■  %  *  (*  -  s)'V?e  •  «25> 

Next,  using  relativistic  kinematics,  we  obtain 

Ws.-  •  «26' 

Using  equation  (26)  in  equations  (23)  to  (25),  then  to  lowest  order  in  it 
and  ic,  we  obtain 

+  (It  -  .  (27) 

Factoring  this  expression  of  energy  conservation  explicitly  into  the  brems- 

strahlung  probability  w?  *  ($,1?) ,  we  define  the  quantity  v?  ±  by 

^a'^p  ^a'*^p 


(28) 


Substituting  equation  (28)  in  equation  (22)  gives 


*  *>£$(<*£  -  **va  +  (£  -  ^)*vg)  . 


Equation  (29)  expresses  the  recoil  force  on  a  particle  of  species  a  due  to 
induced  bremsstrahlung  as  an  integral  over  particle  momentum  j$g,  photon  wave 
vector  and  momentum  transfer  t.  The  integrand  involves  the  particle  dis¬ 
tribution  f|^  of  the  scattering  particle,  the  bremsstrahlung  transition  rate 
v$  ♦  ($,£),  an<J  the  photon  distribution  N-?. 

The  delta  function  in  equation  (29)  expresses  energy  conservation  in  the 
bremsstrahlung  process  in  the  soft  photon  and  small  momentum  transfer  approx¬ 
imation  for  relativistic  particles.  Equation  (29)  agrees  with  Tsytovich’s 
(4.1 exactly.  It  differs  in  sign  from  equation  (3)  of  Akopyan  and 
Tsytovich,2  which  apparently  has  a  typographical  error. 

o 

The  photon  density  N£  in  equation  (29)  can  be  expressed  in  terms  of  the 
associated  field  as  follows.  The  power  density  delivered  to  the  electromag¬ 
netic  field  is  given  by 


du  +  ± 

5 t  ■  -3  E 


V.  N.  Tsytovich,  Bremsstrahlung  of  a  Relativistic  Plasma,  Tr.  Fiz.  Inst. 
Akad.  Nauk  SSSR,  66_  (1973),  191  [Proc  (TRUDY)  of  P.  N.  Lebedev  Physics  Inst., 
66_  ( 1975) ,  199]. 

2A.  V.  Akopyan  and  V.  N.  Tsytovich,  Bremsstrahlung  in  a  Nonequilibrium 
Plasma,  Fiz.  Plasmy,  1_  (1975),  673  ( Sov .  J .  Plasma  Phys.,  1_  (1975),  371]. 
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Using  Maxwell's  aquation 


(31) 

in  aquation  (30)  produces 

it  -  -  -  ll)-£  • 

Equivalently, 

^  -  $.(£xfi)  -  Mxt  +  b  .  (33) 

at  v  J  at 


Substituting  Maxwell's  equation 


$x|  - 


(34) 


in  equation  (33)  produces 


du 

dt 


(35) 


Equation  (35)  is  the  energy  conservation  theoren 


where  I  is  the  Poynting  vector 


8 


t  =  8x8  , 


u  is  the  energy  density  in  the  electromagnetic  field,  and 


3U  S  a£  ±  3$ 

—  =  H*  —  +  £•  —  , 

at  at  at 


The  Fourier  representation  of  the  electric  field  is 


i  -  /  d3£  du> 


&  m  j  d3£  du  Sfce1  (***-<**) 


1 


Also,  by  the  plasma  constitutive  relations 


°ki  ■  eijdt,u)Ekj 


Furthermore,  taking  the  Fourier  transform  of  equation  (34)  produces 


Also, 


t  -  i-  -  eQc2i  . 
^o  ° 


A  A  r«*  .  rm  44/./,  AAA  A  ^  ‘  '  .  *  '  •  * 


V-'.  • 


Next,  using  equations  (39)  to  (43)  and  the  appropriate  Fourier  repre¬ 
sentations  in  equation  (38)  gives 


♦  £oc2  /  a*  *>  a*-  a«*  (*^0i  . 

O  J  m  '  '  (I) 


By  an  ordinary  vector  identity, 


+ 

k*k'E)c*E)ct  -  k'Bjjtk'^B^  . 


Substituting  equation  (45)  in  equation  (44)  gives 


U  -  -i  /  dk  dk 


'  Iw'Bjtieijtit'ju'  >B*. 


i  [  (ic-rtc '  )•?-(«+<»)'  )t] 

®  t 


where  dk  =  d3$  du>.  Symmetrizing  equation  (46)  in  k  and  k'  produces 


Wi 


_2>r  AVI  4 


*1 

61 


■jr.K  .v.y . 


if  -  -  i  J  *  *  Eki“'£ij,1','»,>Ei('j 


si[  (ic+lc • )  •r-(w+u>'  )t] 


The  total  field  involved  in  the  bremsstrahlung  process,  Ek,  is  given 
elsewhere, 5  namely, 


4  =  lg<°>  +  , 


where  ^  is  the  lowest  order  stochastic  bremsstrahlung  radiation  field 
and  is  the  regular  nonradiative  component.  The  bremsstrahlung  field 
written  in  terms  of  its  polarization  vector  is  given  by 


eg<0)  -«10>  - 


where 


«•«*  -  1  • 


SH.  e,  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi- 
librium  Relativistic  Beam-Plasma  System,  Harry  Diamond  Laboratories,  HDL-TR- 
1996  (1983), 


*  v  yv<vv  •: 


■.  > .  • . 
.A'  •  • 


The  stochastic  properties  of  the  brensstrahlung  field  are  approximated  to  the 
needed  order  by 


<Eki0)>  =  0 


(51  ) 


and 


<eki0>t»,(?l>  •  ekiekj  l^10’!2  «(*  +  kl)  152) 

(these  are  eq  (49)  and  (50)  of  a  previous  work5).  Therefore,  substituting 
equation  (48)  in  equation  (47),  keeping  only  the  lowest  order  bremsstrahlung 
field,  integrating  equation  (47)  over  time,  and  using  equation  (51),  we  find 
that  the  ensemble  average  energy  density  in  the  field  is  given  by 

<u>  *  <u>°  +  <u>R  ,  (53) 


where 


<“>°  -  /  at>  -  J  <*  [e1J(t,u>l«<E°<0>E°i°>> 

J  ei  [(£+&'  )»r-(u)+ta'  )t] 


(54) 


eoc2<M£j0)*-.|£(0>> 


D 

is  the  bremsstrahlung  part,  and  <u>  is  the  regular  nonradiative  part.  Only 
<u>°  is  needed  for  the  present  work. 


5ff.  E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi- 
librium  Relativistic  Beam-Plasma  System,  Barry  Diamond  Laboratories,  HDL-TR- 
1996  (1983). 


.  WV.  V.V. 


Substituting  equation  (52)  in  equation  (54)  and  using  the  definition  of 
the  four-dimensional  delta  function. 


6(k)  -  63(*)6(w)  =  (2it)“4  /  d dt  , 


produces 


<u>°  “  \  /  *  dk’  £ekieij  (it»  + 

K,0,l  **(*  *  *■)  • 


In  equation  (56)  we  define 


I(k)  =  /  dk*  «3(Jc  +  £')l«kieij  (*»“)«•&  +  •ki«'eij(ic,#«,)eicj] 


Changing  the  variable  u*  to  ui"  =  w  +  u>‘  and  integrating  over  k',  equation  (57) 


becomes 


I(k)  *  /  — -  6(u")[e£ieij(k,  M)ue£j  +  (a)"  -  a>)e^ieij  (-k,  a)"  -  u>]e£j] 


.  (58) 


Expanding  the  second  term  in  the  integrand  of  equation  (58)  in  a  Taylor 
series  about  a>"  =  0,  then 

100  *  /  “ £  5(“")  ”  wekieij(”^*  ”“)ekj 

'  (59) 

+  «"  -~r  [(“"  ■  w)ekieij(~^' u"  -  • 

By  the  reality  of  the  fields,  it  follows  from  equation  (41 )  that 

-uj)  *  e£j(fc,  u>)  .  (60) 


Substituting  equation  (60)  in  equation  (59)  and  simplifying  the  second  term 
causes  equation  (59)  to  become 


I(k)  =  2i  /  Im  e^*eij (it,  w)e^j 


+  [“ekieij(-*'-a))ekj]  ' 


The  dielectric  constant  ea()t,  u> J  for  mode  a  is  defined  in  terms  of  the 
dielectric  permittivity  tensor  w]  and  the  unit  electric  polarization 

vectors  by 


e°(It,u)  =  egeijAeOegj  +  eQ  (*«*g)(**eg*) 


(see  Other  Works  by  Tsytovich,  Selected  Bibliography). 


Noting  that  the  second  term  of  equation  (62)  is  real,  and  substituting 
equation  (62)  in  equation  (61),  we  obtain 


I(k)  -  2iu>  Im  e°(k,  u)  /  ^  6(w") 


+  3^  [uekieij(-^'“w)ekj]  • 


If  mode  decay  or  growth  are  ignorable,  then  the  mode  dielectric  constant 
c°(k,  (d)  is  real  and  the  first  term  of  equation  (63)  may  be  dropped.  In  that 
case  equation  (63)  becomes 


1<k)  *  [wekieij(-^-w)ekj]  • 


If  we  next  integrate  the  second  term  in  equation  (56)  and  substitute 
equations  (50),  (57),  and  (64),  equation  (56)  becomes 


<u>°  -  /  d*  U 

♦  -*•«**•«)] I l*(0,l 2  . 


•  V  V  V  V  v  -  •  -  .  • .  -  - 


t  »*»  4'  «“*  « .  ♦* 


•••  v  •  •  N-  N* 


Changing  the  variable  of  integration  from  k  to  -k  in  the  first  term  onlyr 
using  the  reality  property  of  the  fields  to  replace  e^E^®^  bY  e£  * 

substituting  equation  (62),  and  simplifying,  causes  equation  (65)  to  become 

o  T3  /  a,,-*-  .  c2  +a‘) 

<u>  *  -  /  dk  j^-gu  (we  (k,  w)  -  e0  —  **ekk*ek  / 

(66) 

Simplifying  equation  (66)  results  in 


<U>°  -  if  dk  [«<><*,«)  +  »  *  E„  c2]|e£I°>|2  . 


(67) 


The  zeroth-order  dispersion  relation  for  the  bremsstrahlung  wave  is  given  by 
the  poles  of  the  photonic  Green's  function  G^j  in  equation  (22)  of  a  previous 
paper5  or,  equivalently,  the  zeros  of  the  determinant  of  G^j  in  equation  (18) 
of  that  paper.5  Thus  the  bremsstrahlung  field  ^  must  satisfy 


- 3 - -  (k.k.  -  k26  .)  +  e.. 

u  (ui  +  i 6 ) 2  13  13  13 


L  o 


E?j0)  -  0 


(68) 


Taking  the  inner  product  of  equation  (68)  with  e^®  and  using  equations  (49) 
(62)  and  y”1  =  eQc2  produces 

(*?*•«?•< t  -  *2)  *  4  (*•*?)(*•*?*)  *  0 


e_c^ 


u)2 


(69) 


or 


k2c2  _ 


W2  eo 


(70) 


SH.  E.  Brandt,  The  Total  Field  in  Collective  Bremsstrahlung  in  a  Nonequi¬ 
librium  Relativistic  Beam-Plasma  System,  Harry  Diamond  Laboratories  HDL-TR- 
1996  (1983). 


Substituting  equation  (70)  in  equation  (67)  results  in 


<u>°  =  ^  J  dk  (2e0(*,  u>)  +  ui  (*,u>))|e£(0)|2  .  (71) 


Equivalently,  equation  ( 71 )  becomes 

<u>°  =  1  ;  d3)t  du>^  (e°(ta>)a>2]|E^(0)|2  .  (72) 

One  can  also  express  the  energy  density  in  the  field  in  terms  of  the  photon 
number  density  as  follows: 

<u>°  =  J  u£n ®  ,  (73) 

(2tt)3  ^  * 


where  uj£  is  the  frequency  of  the  mode  o  as  determined  by  the  dispersion  rela¬ 
tion.  Equivalently,  equation  (73)  may  be  written 

«->°  -  /  *?«(•  -  *£)  • 

Comparing  equations  (72)  and  (74),  we  find 

*ka(“  “  <*£)  *  (e°(^,w)u2)|B^'0)|2  .  (75) 

Equation  (29)  may  be  equivalently  written 
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*3g  do,  d3;  _  JH,  )v{  „  {t  t)  j 

“  j  (2it)9  L  j  1  pe  P&J  Pa'PfT  '  * 


X  6(a)  -  K»va  +  (<  -  1c)*Vg)6(a)  -  u£)  , 


where  the  integral  of  a  delta  function  over  o>  has  been  inserted.  Therefore, 
after  substituting  equation  (75)  in  equation  (76),  we  finally  obtain 


**  ’  -4"3  /  i_  |_  (u2eo,iE,u))i;(j  - 


*  vi<,.e>e(£-*)Kl0,rs(“  -  *•*.  *  (*  -  *Hs) 


Equation  (77)  expresses  the  general  form  for  the  force  cm  a  particle  due  to 
induced  bremsstrahlung.  The  dispersive  properties  of  the  beam-plasma  system 
enter  explicitly  in  equation  (77)  through  the  dielectric  constant  e°(lt,  a))  and 
implicitly  through  the  transition  probability. 


5.  CONCLUSION 


The  general  form  has  been  derived  for  the  collective  bremsstrahlung  recoil 
force  in  a  nonequilibrium  relativistic  beam-plasma  system,  namely. 


-  J  d3pP  a3|t  )v«  .  (tt) 

a  *  (2x)9  v  J  V  Pa'P3 


x  n|6(u£  -  t-ta  ♦  (*  -  *)•♦*)  . 


Here  is  the  force  on  a  particle  of  species  a  and  velocity  va  due  to  its 
participation  in  a  bremsstrahlung  process  in  which  it  interacts  with  a  parti¬ 
cle  of  species  3  and  velocity  $g,  resulting  in  a  momentum  transfer  t  -  it  to 
that  particle.  A  bremsstrahlung  photon  of  wave  vector  $.  and  frequency  <*£  in 

mode  o  is  emitted.  The  functions  f±  and  N$  are  the  particle  and  photon 

o  ® 

distributions,  respectively,  and  V+  is  the  bremsstrahlung  transition  rate 

with  energy  conservation  already  factored  out.  The  photon  distribution  func¬ 
tion  Nj?  can  be  expressed  in  terms  of  the  bremsstrahlung  field  by  equation  (75) 
to  obtain  an  equivalent  form  given  by  equation  (77).  Equation  (78)  was  used 
by  Akopyan  and  Tsytovich1'2  to  obtain  an  expression  for  the  collective  brems¬ 
strahlung  transition  rate  by  means  of  direct  comparison  with  another  expres¬ 
sion  for  the  collective  bremsstrahlung  recoil  force,  which  was  determined  from 
the  equation  of  motion  for  a  dynamically  polarized  test  particle  undergoing 
br ems  s  t rah lung . 

The  present  calculation,  together  with  previous  work  by  the  author  (see 
Selected  Bibliography),  is  important  for  ongoing  work  in  calculating  collec¬ 
tive  radiation  processes  and  conditions  for  the  occurrence  of  radiative  in¬ 
stability  in  relativistic  beam-plasma  systems. 
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